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ABSTRACT 

Lyman-Q- (Lya) emission of neutral hydrogen (/II 215.67 A) is the main contributor to the ultraviolet flux of low-mass stars such 
as M dwarfs. It is also the main light source used in studies of the evaporating upper atmospheres of transiting extrasolar planets 
with ultraviolet transmission spectroscopy. However, there are very few observations of the Lyor emissions of quiet M dwarfs, and 
none exist for those hosting exoplanets. Here, we present Lya observations of the hot-neptune host star GJ 436 with the Hubble 
Space Telescope Imaging Spectrograph (//5r/STIS). We detect bright emission in the first resolved and high quality spectrum of a 
quiet M dwarf at Lya. Using an energy diagram for exoplanets and an A^-body particle simulation, this detection enables the possible 
exospheric signature of the hot neptune to be estimated as a ~ 11% absorption in the Lya stellar emission, for a typical mass-loss 
rate of lO"' g s"' . The atmosphere of the planet GJ 436b is found to be stable to evaporation, and should be readily observable with 
HST. We also derive a correlation between X-ray and Lya emissions for M dwarfs. This correlation will be useful for predicting the 
evaporation signatures of planets transiting other quiet M dwarfs. 
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1. Introduction 

Planetary transits provide golden opportunities to probe the at- 
mospheres of extrasolar planets. The atmospheric signals de- 
tected by this method in the visible (e.g., Charbonneau et al. 
2002) are usually tenuous, while the infrared detections of atmo- 
spheric signatures remain disputed (e.g., Gibson, Pont & Aigrain 
2010). In the ultraviolet (UV), planetary transits of the hot gi- 
ant planet HD 209458b are detected from conspicuous spectro- 
scopic signatures, which are seen as absorption in the emission 
lines originating from the stellar chromosphere and transition re- 
gion between the chromosphere and the corona (see Ehrenreich 
2010 for a review). The strongest signature is a (15 ± 4)% ab- 
sorption in the resolved stellar Lyman-or (Lyor) emission line 
of neutral atomic hydrogen (H i) detected in medium-resolution 
spectra taken with the Space Telescope Imaging Spectrograph 
(STIS) on the Hubble Space Telescope (HST) (Vidal-Madjar et 
al. 2003; see also Ehrenreich et al. 2008). Linsky et al. (2010) 
present new observations with the Cosmic Origins Spectrograph 
(COS) on HST. They report a (7.8 ± 1.3)% absorption in the Cn 
line, confirming the previous result of Vidal-Madjar et al. (2004). 

These distinctive UV signatures, compared to those for the 
visible transit of the whole planet (1.6%; Charbonneau et al. 
2000; Henry et al. 2000), require the presence of an extended 
H I upper atmosphere, or exosphere, to the planet. This envelope 
must be evaporating because (i) it must fill the planetary Roche 
lobe to account for the size of the absorption, and (ii) the H i 
atoms must be accelerated by the stellar radiation pressure be- 
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yond the escape velocity of the planet to account for the absorp- 
tion in the wings of the Lyor line (Vidal-Madjar et al. 2003). The 
presence of elements heavier than hydrogen (C, O, Si) at high 
altitudes requires an hydrodynamic escape (or blow-off) of the 
upper atmosphere, with the flow of escaping H i carrying heavier 
species up to the Roche limit (Vidal-Madjar et al. 2004; Linsky 
et al. 2010), where they are swept away by radiation pressure 
(Lecavelier des Etangs, Vidal-Madjar & Desert 2008). 

New cases of atmospheric escape have been recently re- 
ported for two hot jupiters. Lecavelier des Etangs et al. (2010) 
used HST I ACS to detect the Hi exosphere of HD 189733b at 
Lya. Fossati et al. (2010) observed the highly irradiated planet 
host star WASP- 12 with HST/COS. They interpret the tentative 
detection of extra absorption in the metal lines of the star WASP- 
12 during the transit as signs of mass loss. Thus, atmospheric 
evaporation could be a common phenonemon among close-in 
planets. 

These results have inspired a comprehensive modelling ef- 
fort (e.g., Garcia-Munoz 2007) suggesting that the upper atmo- 
spheres of close-in giant planets must be heated by the stellar 
X-rays, and both extreme and far UV (Cecchi-Pestellini et al. 
2009) to about 10000 K. This temperature has been inferred in 
the thermosphere of HD 209458b from the detection of excited 
hydrogen atoms (Ballester, Sing & Herbert 2007). The typical 
atmospheric mass-loss rate m derived from models is between 
10'" and 10" g s"'. In this case, evaporation should not severely 
impact the stability of giant planets' atmospheres. On the other 
hand, it might be more important for lower-mass planets such as 
hot neptunes and super-earths, if their mass-loss rates are of the 
same order as for more massive planets. 

In this article, we evaluate whether mass loss from a hot 
neptune such as GJ 436b could be detected with current instru- 
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Table 1. Properties of GJ 436 atid its planet 



Parameter 




Value 


Ref. 


The planet 


Mass M,, (Me) 




23.17 ±0.79 


(1) 


Radius R,, (Re) 




4 79+0.09 

^•^^-0.10 


(1) 


Mean density p,, (g cm"') 






(1) 


Gravity gp (m s"^) 




12.8 ±0.8 


(1) 


Semi-major axis a,, (au) 


0.02872 ± 0.00027 


(1) 


Brightness temperature Tb 


(K) 


717 ±35 


(3) 


The star 


Distance dj, (pc) 




10.23 ±0.24 


(4) 


X/EUV luminosity 








log(Lx/Euv [erg s"']) 




26.85+q;8^ 


(2) 


Notes. (1) Torres 2007; (2) Hunsch et al. 1999, Sanz-Forcada et 
al. 2010, Poppenhaeger et al. 2010; (3) Demory et al. 2007; (4) 
Hipparcos catalogue, Perryman et al. 1997. 



mentation, and what its UV transit signature could be. GJ 436b 
is the first transiting Neptune-mass planet discovered (Butler 
et al. 2004). It orbits a close-by (10.2 pc) M2.5v dwarf with 
V — 10.68, and triggered unprecedented interest when Gillon et 
al. (2007a) detected a photometric transit of the planet, inferring 
a true mass of 23.17 M^jand a radius of 4.22 R{5, which are simi- 
lar to those of either Uranus or Neptune. (The planet properties, 
taken from ToiTes (2007), are summarized in Table 1 along with 
some stellar parameters.) However, GJ 436b differs significantly 
from the Solar System ice giants because of its moderately ec- 
centric orbit with a semi-major axis of 0.029 astronomical unit 
(au). At this distance from its 0.026-Lq parent star, the planet 
indeed receives ~ 30 000 times more flux than Neptune receives 
from the Sun, and has a blackbody biightness temperature of 
Tg - in ± 35 K (Demory et al. 2007). Using the Spitzer Space 
Telescope, Gillon et al. (2007b) inferred a planet mean density 
of ~ 1.7 g cm"-', too weak to be that of a body composed ex- 
clusively of water ice, so that GJ 436b must possess a hydro- 
gen/helium envelope that would account for ~ 10% to 20% of 
the planet's mass (Fortney et al. 2007; Figueira et al. 2009). The 
detection of this envelope with transmission spectroscopy could 
validate the internal structure models predicting its presence. 

To assess this possiblity, it is necessary to derived both the 
stellar Lya emission line profile and brightness. We present in 
Sect. |2]the observations that allowed us to calculate the stellar 
emission as seen from Earth. We had to correct for the inter- 
stellar medium (ISM) absorption before estimating the intrin- 
sic stellar Lya brightness: this analysis is presented in Sect. [3] 
In Sect. |4] we estimate a range of plausible mass-loss rates in 
GJ 436b's upper atmosphere. This range serves as an input to 
an atmospheric escape model used to estimate the transit signa- 
ture of the evaporating hydrogen envelope of the planet in the 
UV (Sect. [5]|. To our knowledge, these observations are the first 
resolved accurate measurement of a quiet M dwarf's Lya emis- 
sion. Hence, we discuss in Sect.|6]how they compare to archival 
Lya observations of active M dwarfs. 

2. Observations and data reduction 

We obtained HST time (GO/DD#11817) to precisely estimate 
the Lya emission of the M dwarf GJ 436. The data were 
recorded on 2010 January 6 with the STIS instrument (Woodgate 
et al. 1997). The observation consists of one 1762-s exposure 
on the Far Ultraviolet Multi-Anode MicroChannel Array de- 
tector (FUV-MAMA). The STIS/FUV-MAMA is a solar-blind 
cesium iodide (Csl) detector with a 25" x 25" field of view 
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Fig. 1. Two-dimensional spectral image of GJ 436. The light is 
dispersed along the x-axis (roughly). The stellar Lya emission 
is seen on both sides of the geocoronal emission, which impacts 
the detector across the whole y-axis. 



(~ ^.'025 pixeL'), operating from 1 150 to 1 700 A. The light 
was diffracted using the long slit of size 52" x U! 1 and the first- 
order grating G140M. This grating was tilted to ensure a central 
wavelength of 1 222 A. In this configuration, our data have a 
spectral coverage of 1 194-1 249 A with a medium resolution of 
~ 10000 (30 km s"') and a throughput of between ~ 1.2% and 
2% (see Proffitt et al. 2010). The data were recorded in time-tag 
mode, and treated with the version 2.27 of CALSTIS, the STIS 
data pipeline. CALSTIS performs all the standard data reduction 
tasks (such as bias and dark subtraction, flat fielding, and wave- 
length and flux calibrations; see Dressel et al. 2007 for more 
details). For spectroscopy, the pipeline end-products for a single 
exposure are the two-dimensional spectral image and the one- 
dimensional spectrum. 

Figure[T]shows the geometrically rectified, wavelength-, and 
flux-calibrated 2-D spectral image (X2D). The spectral disper- 
sion is along the x-axis. In this image, the stellar Lya emission 
appears on both sides of the prominent geocoronal emission, 
which is imaged across the whole detector y-axis. This terres- 
trial air glow seriously degrades the quality of spectroscopic ob- 
servations at Lya; however its rather naiTow extent along the 
dispersion (x-)axis and its large spatial extent (y-axis) enable us 
to perform an efficient correction. This task is handled by the 
pipeline subroutine called BACKCORR, which performs the 1-D 
spectral background subtraction. The background is determined 
along a five-pixel-wide region located +300 pixels away from 
the centre of the spectral extraction region. The 1-D extracted 
background is then fitted by a «-th degree polynomial, which re- 
places the background except in the regions surrounding the Lya 
line (Dressel et al. 2007). The resulting calibrated 1-D spectrum 
of GJ 436 (XID) is shown in Fig.|2] 

3. Lyman-a emission of GJ 436 

3. 1 . Observed flux at Earth 

Figure |3] shows a zoomed region of the recorded spectrum 
around the Lya line. The error bar for each spectral bin is es- 
timated by CALSTIS, which basically propagates the statistical 
noise ■\f{Ng)/g in each pixel of the raw data, where is the 
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Fig. 2. The reduced and calibrated STIS/G140M spectrum of 
GJ 436 in January 2010. 



2000; Hebrard & Moos 2003; Linsky et al. 2006) that this value 
of the D/H ratio is constant within the Local Bubble (< 100 pc; 
see Wood et al. 2005 and references therein). 

The ISM column density value is chosen within the range 
considered by Wood et al. (2005), 17.6 < log NiH i) < 18.2. The 
value is adjusted to 10"* cm"^ using the D i line at 1 215.25 A. 
This agrees with the ISM column density of logA^(Hi) = 17.82 
measured by Wood et al. (2005) toward the star HD 4345, in 
a similar direction (11.1° from GJ 436) and at a distance of 
21.7 pc. A column density logA^(Hi) » 18 would be unreal- 
istic for a target as close as GJ 436 (10.2 pc) and yield a stronger 
D I absorption than observed. (The presence of this feature can 
be guessed in the blue part of the observed Hi Lya emission 
in Fig. |3]) In any case, the fit to the observed Lya profile does 
not depend much on the assumed ISM column density between 
A^(Hi) ~ 10'^ and 10'^ cm"^ because at these column densities, 
the H I absorption line is saturated. The result is more dependent 
on the assumed profile of the wings extrapolated to the core of 
the emission line. 



number of data events (counts) and g is the gain (equal to 1 pho- 
ton count"' for MAMA observations). 

The total flux of the line seen by HST, integrated from 
1 214.40 to 1 215.58 A and from 1 215.75 to 1 217.00 A in or- 
der to avoid any residuals from the air glow subtraction at the 
centre of the Une, is (5.490 ± 0.184) x lO"'"* erg s"' cm"-. We 
note that, as seen from the Earth, the Lya flux of GJ 436 is higher 
than that of HD 209458b: the stellar emission from the latter star 
is also reproduced in Fig. |3]for comparison. 

The observed shape of the line is a characteristic double peak 
resulting from the addition of several profiles, as shown in Fig. [4] 
As for the Sun (Woods et al. 2005), the core of the intrinsic stel- 
lar Lya line is an emission feature originating in the transition re- 
gion, where the temperature increases quickly between the chro- 
mosphere and the corona. The wings of the Lyor line originate 
in a deeper chromospheric region. The emission can be either 
single-peaked or double-peaked, with a central dip due to the 
high opacity of the abundant hydrogen atoms. As seen from the 
Solar System, the central part of the line is strongly absorbed 
by the ISM from 1 215.58 to 1 215.75 A. The ISM absoiption 
can indeed be Doppler-shifted relative to the line centre by sev- 
eral kilometers per second. The ISM absorption is the combina- 
tion of absorption by neutral hydrogen (H i) and deuterium (D i). 
Deuterium produces narrow absorption at 1 215.3 A, which is 
blue-shifted by ~0.3 A from the Hi absorption line. The ob- 
served Lya profile is therefore caused by the stellar emission line 
being absorbed by the ISM and subsequently convolved with the 
instrumental line spread function (Fig.|4]). 

3.2. Reconstructed stellar intrinsic emission 

To model the hydrogen atom dynamics under the stellar Lya ra- 
diation pressure, we need to estimate the stellar emission line as 
seen by hydrogen atoms escaping from the planet. To calculate 
the resulting Lya transit light curve, we also need to estimate the 
absorption by the escaping atoms over the stellar profile as seen 
from the Earth (Sect. |5]l. To produce these two line profiles, 
we fitted the observed Lya profile with a stellar emission line 
and an ISM absorption model, following the method proposed 
by Wood et al. (2005). In the following, we assume a D/H ratio 
of 1.5 X 10"^. There is widespread consensus (e.g., Ferlet et al. 



As pointed out by Wood et al. (2005), the reconstruction of 
the intrinsic Lya emission of a star can be affected by absorp- 
tion from both the heliosphere and the astrosphere, which we did 
not model for GJ 436. Owing to the moderate resolution of the 
present observations and our lack of knowledge about the ISM 
absorption, from either D i or heavier elements (Fe ii or Mg ii), 
this absorption remain unconstrained. However, GJ 436 is lo- 
cated 87.7° from the upwind direction of the ISM flow seen 
by the Sun, hence based on Fig. 13 in Wood et al. (2005) it is 
unlikely that heliospheric absorption is detectable. No extra H i 
absorption is required besides the ISM component to achieve a 
good fit of the Lya line. In particular, we did not consider an as- 
trospheric absorption component. Wood et al. (2005) found that 
this absorption was necessary in some case to obtain satisfactory 
fits to the Lya lines. This does not seem to be the case for GJ 436 
given the current data resolution. 

We performed the model fitting using several intrinsic line 
profiles, including a single-peak Gaussian and a double-peaked 
line modelled by a double Gaussian. We found that the double- 
peaked Gaussian reproduces far more accurately the observed 
wings in the line profile. The model intrinsic profile is composed 
of two Gaussians with the same full width at half maximum 
(FWHM = AA) and line centres separated by the same value, 
AA, which is taken as a free parameter. The line total flux is the 
second free parameter. The ISM and the stellar radial veloci- 
ties are also free parameters. The model is convolved with the 
G140M instrumental line spread function, and compared to the 
observations. We found a satisfactory fit with AA - 0.41 A, an 
ISM radial velocity of -4 km s ' in the stellar reference frame, 
and a stellar radial velocity corresponding to a Doppler shift of 
0.003 A (~ 0.7 km s"'). The resulting is 34.5 for 39 degrees 
of freedom (see the data and fit in Fig.|4]l. 

The total flux in the reconstructed Lya line is (2.7 ± 0.7) x 
10"'^ erg s"' cm"^. The error bars were estimated by exploring 
the parameter space for the fit to the observed H i line wings, and 
considering possible single-peaked profiles or the possibility of 
two narrow peaks with a deep self-absorption in the central part 
of the emission line. These last two cases produce extreme Lya 
emission profiles within the ISM absorbed wavelength range; the 
resulting error bars should be considered as conservative. 



Ehrenreich, Lecavelier des Etangs & Delfosse: HST/STIS Lyo- observations of GJ 436 

n — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




1214 



1215 . 1216 

Wavelength (A) 



1217 



Fig. 3. Zoom on the Lya line of GJ 436 measured with 
STIS/G140M with error bars from the STIS pipeline (blue 
curve). The STIS spectmm of HD 209458 obtained by Vidal- 
Madjar et al. (2003) is shown for comparison (grey curve). As 
seen from Earth, the nearby M dwarf clearly shows a more in- 
tense stellar emission and a narrower ISM absorption. The lo- 
cation of the air glow (subtracted here) is indicated by green 
hatches. 



4. Mass-loss rate of GJ 436b 

The mass-loss rate ih of GJ 436b is a free parameter in our mod- 
elling of the evaporation signature described in the next section. 
Nevertheless, estimating its value is useful to constraining all the 
parameter space. To achieve this, we use the method presented 
by Lecavelier des Etangs (2007), which consists of locating the 
planet in an "energy diagram". The energy diagram of exoplan- 
ets measures the whole gravitational potential energy of a planet, 
as a function of high-energy irradiations it receives. This dia- 
gram was updated by Davis & Wheatley (2009), and is being 
completed with a larger sample of transiting planets (Ehrenreich 
& Desert 201 1), including GJ 436. 

The resulting diagram is presented in Fig. |5] For GJ 436, we 
calculate, following Lecavelier des Etangs (2007), the potential 
energy per mass unit of the planet 



dE' 



dm 



^ = -GMpIRp + brides = -2.05 X 10'2 erg g- 



(1) 



where G is the gravitational constant and ^tides is the modifica- 
tion of the planet gravitational potential by the stellar tides (see 
Appendix B of Lecavelier des Etangs 2007). 

The energy received by unit of time by GJ 436b can be ex- 
pressed as 



dE 



X/EUV 



df 



-KpGp i^x/EUV 



2.19 X 10^^?7 erg s 



(2) 



The 'X/EUV in Lx/euv means that the luminosity should be 
integrated from ~ 1 A to 912 A, following Cecchi-Pestellini et 
al. (2009) {^However, observations do not usually cover a suf- 
ficiently wide bandpass to yield a precise value of L ~' 



'X/EUV 



The 



' In this paper, we use the following conventions for naming spec- 
tral regions: X-rays (X) in 1-10 A, soft X-rays (XUV) in 10-30 A, 
extreme ultraviolet (EUV) in 30-912 A, and far ultraviolet (FUV) in 
912-1216 A. 
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Fig. 4. Plot of the theoretical profile of GJ 436's Lya line. The 
black thin line shows the intrinsic stellar emission line as seen by 
hydrogen atoms escaping the planetary atmosphere. The black 
thick line shows the resulting profile after absorption by the ISM 
hydrogen and deuterium. The line profile convolved with the 
HST G140M instrumental line spread function (red line) is com- 
pared to the observations (blue histogram), yielding a good fit 
with a.x^ of 34.5 for 39 degrees of freedom. The location of the 
(subtracted) air glow is indicated by green hatches. 



Rosat luminosity (Hiinsch et al. 1999) given in Table[T]is calcu- 
lated from 0.1 to 2.4 keV, i.e., from 5 to 124 A. It is, therefore, a 
lower limit to the actual X/EUV luminosity. 

The factor tj is the heating efficiency in the exoplanet ther- 
mosphere. Lecavelier des Etangs (2007) considered the extreme 
case of // = 1 where all the stellar flux is used to escape the atmo- 
sphere. Considering (without quantifying) energetic losses due 
to thermal emission by atmospheric hydrogen, Tian et al. (2005) 
modelled the atmospheric escape process of the hot Jupiter 
HD 209458b by assuming that rj - 0.15. This value was initially 
chosen by Watson, Donahue & Walker (1981) in their pioneer- 
ing study of Earth's atmospheric escape. In the following, we 
consider both of these values of rj. 

The resulting mass-loss rate is thus 



dE 



X/EUV 



dE' 



(3) 



For // = 1 and 0.15, we find that m = 1.07 x lO"* and 
1.60 X 10^ g s"', respectively. These values are close to the 
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Fig. 5. Energy diagram for 100 transiting exoplanets, updated 
from Lecavelier des Etangs 2007. The energy needed to escape 
a unit of mass of the planet atmosphere is plotted versus the 
X/EUV flux. The dotted lines indicate constant mass-loss rates. 
The blue dot indicates the location of GJ 436b, calculated here 
assuming an heating efliciency of 77 = 0.15. The horizontal error 
bars represent variations in 77 of between 0.01 and 1. 



dM/dt=10° g/s 




Velocity (km s ') 

Fig. 6. The GJ 436 Lya spectrum during the transit of GJ 436b. 
The star spectrum before or after the transit is shown by a thick 
grey line. The spectrum resulting from the absorption by the 
planetary exophere during the transit is shown for various es- 
cape rates. The blue dotted, the red dot-dashed, and black long- 
dashed lines show the resulting spectrum when the escape rate is 
10^, 10^, and 10'° g s"', respectively. With this last escape rate, 
the blue side of the line is absorbed by about 17%, resulting in a 
decrease of about 1 1 % in the total Lyman-a flux. 



canonical escape rate predicted by several evaporation models 
(e.g., Yelle 2004, 2006; Lecavelier des Etangset al. 2004; Garcia- 
Munoz 2007; MuiTay-Clay, Chiang & Murray 2009). 

5. Observable signature of the evaporating 
atmosphere 

To calculate a theoretical H i transit light curve for various atmo- 
spheric escape rates, we modelled the atmospheric gas escaping 
GJ 436b with a numerical simulation including the dynamics of 
the hydrogen atoms. In this A^-body simulation, hydrogen atoms 
are released from GJ 436b's upper atmosphere as particles with 
a random initial velocity corresponding to a 10000-K exobase. 
This temperature corresponds to the temperature expected in the 
upper atmospheres of hot jupiters (e.g., Lecavelier des Etangs 
et al. 2004; Stone & Proga 2009) and was inferred from ob- 
servations (Ballester, Sing & Herbert 2007; Vidal-Madjar et al. 
201 1). In any case, our results do not depend on this assumption 
because atoms are rapidly accelerated by the radiation pressure 
to velocities several times higher than their initial velocities. In 
this model, the gravity of both the star and the planet are taken 
into account. The radiation pressure from the stellar Lya emis- 
sion line is calculated as a function of the radial velocity of the 
atoms. The self-extinction within the cloud is also taken into ac- 
count. The hydrogen atoms are supposed to be ionized by the 
stellar EUV flux (< 912 A). In the simulation, their lifetime 
is calculated as a function of the ionizing EUV flux which is 
taken to be equal to the solar value. The only free parameter is 
the atomic hydrogen escape rate, expressed in grams per second. 
The dynamical model provides a steady state distribution of po- 
sitions and velocities of escaping hydrogen atoms in the cloud 
surrounding GJ 436b. From this information, we calculated the 
corresponding absorption over the stellar line (Fig. |6]) and the 
corresponding transit light curve of the total Lya line (Fig. |7]i. 
We found that escape rates of ~ 10'' to 10'" gs"' can produce 
absorptions of between ~ 3% and 11% in the Lya transit light 
curve when integrating over the whole line. Absorption depths 
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Fig. 7. Plot of the theoretical light curve of the total Lya flux for 
various escape rates m: 10^ g s"' (blue dotted line), 10^ g s"' 
(red dot-dashed line), and 10'" g s ' (black dashed line). The 
thin long-dashed line shows the light curve for the planet without 
additional absorption by the planetary atmosphere. 



this large could be promptly detected with HST Lya observa- 
tions. 

This contrasts with the atmospheric signatures expected in 
the infrared. The atmospheric scale height of this planet is // = 
^^TsKiigp) ~ 232 km, where fi - 2 g moL' is the molar mass 
of a H2 atmosphere and gp is the surface gravity of the planet. 
We can expect atmospheric transit signatures at the level of ~ 
2iAF/F)iH/Ri,) (Winn 2010), where AF/F ^ {RplR^^f ~ 0.7% 
is the transit occultation. This yields a typical extra absorption 
value of ~ lO""' for one atmospheric scale height. This value 
is consistent with the upper limits on the atmospheric absorp- 



6 



Ehrenreich, Lecavelier des Etangs & Delfosse: HST/STIS Lyo- observations of GJ 436 



tion of GJ 436b set by Pont et al. (2009), who have covered two 
near-infrared transits of the planet with the NICMOS camera on 
HST using a grism to sample the 1.1-1.9 /im band. These au- 
thors report observing a flat transmission spectrum at the level 
of ~ 10""*, and in particular no significant signal in the 1.4-yum 
water band (see also Gibson, Pont & Aigrain 2010). 

Using eclipses of the planet by the star observed by Spitzer, 
Stevenson et al. (2010) obtained planet-to-star flux ratios that 
show marginal variations between 3.6 and 24 /im (see also 
Beaulieu et al. 2011). The interpretation of these broad-band 
spectrophotometry data is difficult and relies on statistical ap- 
proaches (Madhusudhan & Seager 201 1). Shabram et al. (201 1) 
modelled the infrared transmission spectrum of GJ 436b. These 
authors conclude that the detailed infrared characterization of 
this planet may be possible using the James Webb Space 
Telescope. Before its launch, the use of UV facilities should be 
considered as a powerful way of characterizing the upper atmo- 
spheres of low-mass exoplanets. 



6. Lya emission of M dwarfs and stellar activity 

Stellar Lyor emission lines are important spectral features in 
the context of exoplanet stellar environment and stellar physics. 
These emission lines are the main contributors to the flux of low- 
mass stars at short wavelengths, from X to FUV. For instance, 
the Lya emission of the Sun contributes to more than 50% of its 
flux below 1 216 A (Woods et al. 1998). The Lya line is used as 
a proxy for determining the temperature and pressure profiles of 
upper stellar atmospheres. 

The Lya emission originate in the transition regions between 
the chromospheres and coronse, while other commonly used ac- 
tivity diagnostics are based on indicators originating from the 
coronse (X-ray flux) or the chromospheres (Mg ii h and k lines). 
In this context, stellar Lya emission provides an almost unique 
access to the structure of the transition region, where the tem- 
perature profile exhibits a deep minimum. Houdebine & Doyle 
(1994) modelled the Lya emission of M dwarfs. According to 
their work, the transition region of these stars is thinner than that 
of solar-like stars. 

Stellar Lya measurements are scarce in the literature. This 
is mainly because it is impossible to observe this emission from 
the ground and that it is difficult to correct measurements for the 
effects of the ISM absorption, both astrospheric and heliospheric 
absorption, and geocoronal emission. Furthermore, correlations 
between Lya emission and other proxies of stellar activity can 
be misleading because they do not originate in the same region. 
This ensures that an accurate estimate of the Lya brightness from 
an activity relation is rather hazardous. In addition, the Lya line 
is extremely variable with time. For the Sun, the integrated Lya 
line flux may change by 37% during one rotation and up to 50% 
over a couple of years (Vidal-Madjar 1975). During the magnetic 
cycle, the extreme (single day) values can vary by more than a 
factor of two (Woods et al. 2000). 

Very few M dwarfs have Lya measurements. Landsman & 
Simon (1993) used the International Ultraviolet Explorer iJUE) 
to report Lya fluxes corrected from the ISM absorption for 12 
M dwarfs. Wood et al. (2005) obtained resolved Lya measure- 
ments for four M dwarfs with //5r/STIS. Three objects are com- 
mon to both samples, leading to a total of 13 measurements. 
This sample consists mostly of active stars. In particular, all the 
HST measurements of Wood et al. (2005) were obtained for M 
dwarfs with magnetic activity levels well above that of GJ 436. 
Therefore, our measurement of GJ 436's Lya flux is invaluable 



for estimating the properties of the transition region in quiet M 
dwarfs. 

The measurements found in the literature are listed in 
Table [2] In Fig. [8] we have plotted the normalized Lya flux, 
^Lyo- = fhya/fhoi ^s a function of the normalized X-ray flux 
Rx - fx/fboi for the lUE and HST measurements from Table |2] 
We have not included the available data for YY Gemini and 
DT Virginis because these stars are close binaries for which a 
definitive identification of each component's flux is difficult. The 
values are listed in Table[2] The observed X-ray fluxes fx are ex- 
tracted from the Rosat all-sky survey (RASS) in the Nexxus 2 
database (Schmitt & Liefke 2004). The bolometric fluxes /boi 
are obtained by combining the IK photometry of Leggett (1992) 
with the bolometric correction BC^ versus I - K relation of 
Leggett et al. (2000). 

As can be seen in Fig. |8] the Lya fluxes extracted from HST 
are significantly lower - for a given Rx - than those from lUE. 
Wood et al. (2005) note that the HST and lUE data are coherent 
for the majority of the objects (besides M dwarfs), except for the 
fainter ones for which the lUE fluxes are systematically higher 
than the HST ones. They attribute this to a bias in the corrections 
for the ISM absorption and geocoronal emission, which are par- 
ticularly difficult at the lower resolution of lUE. 

In the case of faint M dwarfs, we prefer to consider 
the clearly resolved HST measurements, which were obtained 
mainly for magnetically active stars. The values of Ri,ya seem to 
be correlated to Rx for M dwarfs, with log/^Lyi, ~ 0.5 \ogRx - 
2.2, although additional measurements are necessary to confirm 
this trend. The trend implies that Rhyn/Rx > 1 for quiet M 
dwarfs, in contrast to the results for more active ones. Figure |8] 
can thus be used -with caution- to estimate Ri^ya for all M 
dwarfs that have values of Rx between -5 and -3, i.e., cover- 
ing the domain of magnetic activity between 20% of that of the 
quietest M dwarfs in the solar neighbourhood to the very ac- 
tive stars close to the X-ray emission saturation level (see for 
instance Fig. 5 in Delfosse et al. 1998). 

Meanwhile, we stress that the X-ray and Lya fluxes have 
been measured here at different times and, therefore, for dif- 
ferent magnetic states of the stars, which can add dispersion 
to any relationship. For GJ 436, the measured logLx ranges 
from 25.96, as measured with XMM-Newton by Sanz-Forcada 
et al. (2010), to 27.15, as found in the Rosatl?SPC archives by 
Poppenhaeger et al. (2010). The value we have used throughout 
this work comes from the Rosat all-sky survey catalogue of the 
nearby stars (Hiinsch et al. 1999) and is in-between these ex- 
treme values. The dispersion (standard deviation) in these three 
values is cr^^ - 0.77. 

We have retrieved all available X-ray measurements for the 
stars in Table |2] from the Nexxus 2 database and calculated 
their cri^. However, these calculated dispersions only provide a 
rough idea of the real intrinsic X-ray variabilities of the listed M 
dwarfs, because the number of measurements (also reported in 
Table |2]) is not equivalent for all stars. In addition, the measure- 
ments are derived using different instruments, mainly Rosat or 
XMM-Newton (or both). Differences between surveys, e.g. in the 
bandpasses, are not taken into account, thus adding to the disper- 
sion. Finally, the occurrence of a flare during a given observation 
cannot be excluded, leading to a non-gaussian distribution of the 
measurements. Therefore, these values should only be regarded 
as indicative. 
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the Lyo- transit light curve. These figures will allow the existence 
of such an extended upper atmosphere to be tested with a dedi- 
cated HST program. These results also suggest that other quiet M 




dwarfs are brighter at Lya than in X-rays. This prediction opens 
new perspectives for the atmospheric characterization of other 
hot-neptune and even super-earth atmospheres. Interestingly, the 
nature of GJ 1214b, a super-earth, or 'small-neptune' also orbit- 
ing a quiet M dwarf (Charbonneau et al. 2009), could be resolved 
with this technique, provided the star is bright enough at Lya. 
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Fig. 8. Normalized Lya flux Ri^y^ - fiyalfboi as a function 
of normalized X-ray flux Rx - fx/fboi for all the lUE (o) 
and HST (•) Lyo' flux measurements of M dwarfs in the lit- 
erature. Our measurement for GJ 436 data is shown in blue. 
The dotted lines link HST and lUE mesasurements for ob- 
jects that have been measured with both instruments. The plain 
line show the best-fit linear regression to HST measurements, 
log(/?Ly(r) ~ 0.5 log(/?x) " 2.2. The grey dashed line indicates the 
log 7?Lya = log Rx relation. The horizontal error bars in the X-ray 
flux are the dispersions reported in Table |2] 



Table 2. NormaUzed Lya and X-ray fluxes of M dwarfs. 



Star 


Sp.T. 




Rx' 


o-L.^ (dex) 


Nk 




Lyo- 


from HST 








Proxima Cen 


M6v 


-3.88 


-3.58 


0.31"' 


10 


AD Leo 


M4.5v 


-3.50 


-3.12 


0.12"' 


6 


EVLac 


M3.5v 


-3.77 


-2.70 


0.41"' 


8 


AUMic 


Mlv 


-3.51 


-2.97 


0.15"' 


3 


GJ436 


M2.5v 


-4.60 


-4.87 


0.77§ 


3 




Lyo- 


from lUE 








AXMic 


MOv 


-4.22 


-5.29 


0.14"' 


3 


DT Vir 


M2v 


-3.27 


-3.15 


0.41"'"* 


5 


GJ887 


M2v 


-4.05 


-5.13 


0.15"' 


2 


GJ411 


M2v 


-4.44 


-5.05 


0.35* 


4 


AD Leo 


M4.5v 


-3.58 


-3.12 


0.12* 


6 


EVLac 


M3.5v 


-2.74 


-2.70 


0.41* 


8 


YZCMi 


M4.5v 


-3.21 


-3.06 


0.03*-* 


2 


V1216Sgr 


M3v 


-3.34 


-4.13 


1.10*-* 


3 


V998 Ori 


M3.5v 


-2.95 


-3.08 


0.18* 


2 


Proxima Cen 


M6v 


-3.58 


-3.58 


0.31* 


10 



Notes. (*) Normalized X luminosities are Rosat values ex- 
tracted from the Nexxus 2 data base. The dispersion cr^j. is 
the standard deviation in the Nx measurements found in the 
Nexxus 2 database. These measurements are derived using (t) 
Rosat, ($) XMM-Newton, or (t, i) both. (§) The dispersion in 
the case of GJ 436 is obtained as the standard deviation in the 
three Lx values listed in the text. 



7. Conclusion 

We have obtained the first UV glance at a star hosting a tran- 
siting hot neptune. The HST data unambiguously show that the 
early and quiet M dwarf GJ 436 has bright Lya emission. With 
an escape rate of around ~ 10"^ g s"', the hydrogen upper at- 
mosphere of GJ 436b should cause large absorption (~ 11%) in 
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